Mixtures of Bose gases confined in a ring potential 
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The rotational properties of a mixture of two distinguishable Bose gases that are confined in a 
ring potential provide novel physical effects that we demonstrate in this study. Persistent currents 
are shown to be stable for a range of the population imbalance between the two components at low 
angular momentum. At higher values of the angular momentum, even small admixtures of a second 
species of atoms make the persistent currents highly fragile. 
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Introduction. One of the most fascinating phenomena 
associated with superfluidity [T] is the stability of per- 
sistent currents. In some remarkable experiments that 
have been performed recently, Bose-Einstein condensed 
atoms were confined in annular traps 0, [1] > which per- 
sistent currents could be created and observed Q . In an 
earlier experiment, the resistant-free motion of an object 
through a Bose-Einstein condensate below some critical 
velocity, was also observed Q. 

Motivated by these recent advances, in the present 
study we consider a mixture of two (distinguishable) Bose 
gases at zero temperature @, 0] , that are confined to one 
dimension with periodic boundary conditions, i.e. in a 
ring potential, deriving a series of exact and analytic re- 
sults. 

The main issue of our study concerns the rotational 
properties of this system and the stability of persistent 
currents. In higher dimensions it has been argued that 
mixtures of Bose gases do not support persistent cur- 
rents, because there is no energy cost for the system to 
get rid of its circulation (i.e., the line integral of the ve- 
locity field around a closed loop that encircles the ring), 
as long as angular momentum can be transferred between 
the two species 8] . Here, we demonstrate that when the 
total angular momentum per atom varies between zero 
and unity, currents are stable for a certain range of the 
ratio of the populations of the two species. We calculate 
the critical strength of the coupling for a given value of 
this ratio, which we determine analytically and exactly. 
On the other hand, for higher values of the angular mo- 
mentum per atom, persistent currents in one-component 
systems are very fragile, as even small admixtures of a 
second species of atoms destabilize the currents. 

Model. Assuming a ring potential (which corresponds 
to a very tight annular trap along the transverse direction 
@), the Hamiltonian of the system that we study for 
the two components that we label as A and B \s H — 

Haa + Hbb + tjAB Ettf=i '^(^' - ^j )' ^^^''e 
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with k = A,B. 
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2'Kh?'aAB/{MABRS) are the matrix elements for zero- 
energy elastic atom-atom collisions (all assumed to be 
positive), with Mae = MaMb/{Ma + Mb) being the 
reduced mass. Also, R is the radius of the annulus and 
S its cross section, with R ^/S. 

We start from the mean-field approximation, introduc- 
ing the order parameters of the two components (pA and 
(pB', later we also go beyond the mean-field approxima- 
tion, diagonalizing the Hamiltonian H numerically and 
analytically. The resulting (coupled) nonlinear Gross- 
Pitaevskii-like equations are 



5^ 



NkUkk\<l>k\^<l>k + NiUki\<l>i\'^(l)k = /^/=<^/c,(2) 



where / \(j>k\'^d9 — 1. Here ^k are the chemical potentials 
divided by the kinetic energy e — h? /{2MR^), where 
we have assumed for simplicity equal masses for the two 
species. Ma = Mb = M. Also, Uki = Uu/e., with k,l = 
A,B. 

Energetic .stability, dynamic stability and phase sepa- 
ration. Before we turn to the rotational properties, let 
us consider briefly the question of phase separation. In 
homogeneous systems it has been shown that the condi- 
tion for e nerg etic stability of the homogeneous solution 
is [13, [O, m U\b - UaaUbb < 0, and also Uaa > 0, 
Ubb > 0. One may generalize this result for the case 
of a finite system, taking into account the contribu- 
tion of the kinetic energy. The details of this calcula- 
tion will be reported elsewhere. Here we just mention 
that this more general condition is Jab ~ IaaIbb < 
1/4-1- {'^AA + lBB)|'2^ where we have introduced the pa- 
rameters ^k,i = Uk,i^/ NkNi/ [211) for convenience (these 
parameters give the ratio between the typical interaction 
energy and the typical kinetic energy). As one crosses the 
phase boundary, the two clouds develop sinusoidal vari- 
ations in their density, with an amplitude that increases 
continuously from zero. 

The dynamic stability of the system may be ex- 
amined with use of the (two coupled) Bogoliubov-de 
Gennes equations. Again, the details of this cal- 
culation will be reported elsewhere. The dispersion 



Here Mk are the atom 
inh'^akk/iMkRS) and Uab = 



that one obtains from this analysis is lu'^ 

ilAA + IBB ± \/{lAA - IBb) 
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quirement of a real uj implies the same condition as that 
for energetic stability. 

Effect of the periodicity on the dispersion relation. The 
one-dimensional motion that we have assumed in our cal- 
culation, in combination with the periodic boundary con- 
ditions have some important consequences on the disper- 
sion relation, which are also present in the case of a single- 
component gas, as shown by Bloch [13 ]. The matrix el- 
ements that determine the interaction energy do not de- 
pend on the quantum numbers of the angular momentum 
m, and also the center of mass coordinate separates from 
the relative coordinates. As a result, solving the problem 
in the interval < I < 1, where I = (La+Lb)/ {Na+Nb) 
is the angular momentum per particle, then exciting the 
center of mass motion, we may evaluate the spectrum at 
any other interval n < I < n + 1. More specifically, if 
(l>A.o = J2m Cm$m, and (f>B,o = J2m drn^m, arc the order 
parameters for < / < 1, then the order parameters for 
n < I < n + 1 axe given by (j)A,n = Cm^'m-t-n, and 

Denoting the energy per atom for n<l<n+las 
En{l)/N, then E„{l)/N = Eq{Iq)/N + + 2n?o, where 
< /o < 1, and Z = /o + ^- Therefore, En{l)/N - P = 
Eo{Io)/N—Iq, which are both equal to a periodic function 
e{l), i.e., e(lo + n) — e{lo). Thus, we write quite generally 
that 

En{l)/N = P + e{l) = {lo + nf + e{h). (3) 

In other words, the energy of the system for n < I < 
n + 1 consists of an envelope part, i.e., the first term 
on the right, which arises because of the center of mass 
excitation, plus a periodic part e(Z). 

Furthermore, the function e{lo) is symmetric around 
Iq = 1/2 (an example of this symmetry is demonstrated 
below, where it is shown that Eq/N is linear for < I < 
xb = i — XA and xa < I < i)- To see this, let us consider 

the states (j)^ = I]mCm$l-m, and 0f = X^mCim^l-m, 

with an V equal to 1 — I, or I + I' = 1. It turns out that 
the difference in the energy per particle in the states 0^, 
(pg, and is AE/N — I' — I. However, according 

to Eq. (ISD, AE/N = I' - I + e{l') - e{l), and therefore 
e(Z') = e{l), which means that e{lo) is indeed symmetric 
around Iq = 1/2. 

Rotational properties. Since, according to what 
was mentioned above, the dispersion relation is quasi- 
periodic, in order to study the rotational properties of 
the gas, we restrict ourselves to the interval < I < 1. 
We introduce the variables xa — Na/{Na + Nb) and 
Xb = Nb/ {Na + Nb), and assume without loss of gener- 
ality that Xb < Xa, with xa + xb = In what follows 
we also assume equal scattering lengths, and therefore 
Uaa — U BB = Uab = U. The condition of equal scatter- 
ing lengths is not far from reality, with Rubidium atoms 
in different hyperfine states being an example. Interest- 
ingly, in this case there is a series of exact, analytic re- 
sults. If this condition is weakly violated, the deviations 
from these results will be small. 

According to the result mentioned earlier, for Uaa = 



Ubb = Uab the gas is in the homogeneous phase, and 
it is both dynamically, as well as energetically stable. In 
this case, we find that ior < I < xb and xa < I < ^, 
only the states with <I>o and $1 are (macroscopically) 
occupied. The interaction energy of the gas is equal to 
that of the non-rotating system, since the total density 
n{0) = nA{0) + nB{0) is homogeneous. As a result, the 
total energy of the gas varies linearly with I. These are 
exact results within the mean-field approximation. On 
the other hand, for xb < I < xa more states contribute 
to the order parameters, while the dispersion relation 
is not linear in this interval. More specifically, let us 
consider the states of some fixed expectation value of the 
angular momentum I, (f>A,o — co$o + ci^i, and (pB.o ~ 
do$o + with x^lcip +XB\d\^ = I, and also |coP -I- 

|cip = 1, |doP -I- — 1. The above states have a 

maximum value of / equal to unity. Evaluating the total 
energy Eq and minimizing it, it turns out that 

Eq/N ^l + -i{\/2 + (xA\cQ\\c^\ - XB\dQ\\dx\f \ , (4) 

where N — Na + Nb is the total number of atoms and 
7 = NU/{2t:). ¥or Q < I < xb and xa < I < I, 
the last two terms may be set equal to each other, 
which means that Eq/N = I + j/2. Remarkably, any 
other single-particle state cannot lower the energy and 
its occupancy is exactly zero. The occupancies of the 
single-particle states with m = and m — 1 are Cq — 
{xA-l)il-l)/[xAil~2l)], and cj = l{xB-l)/[xAil~2l)]; 
d^ and d^ are given by similar formulae, with xa and xb 
interchanged. The same expressions hold for a mixtures 
of two Bose gases that are confined in harmonic traps 
[l3 |. but in this case the energy is parabolic and not lin- 
ear in I. 

Persistent currents. Let us now examine the question 
of stability of persistent currents. In the case of only one 
component, for 7 > 3/2 , the system supports persistent 
currents at Z = 1 [islll^. As we saw earlier, if one starts 
with XA = ^ and xb = and increases the population 
of the B component, the dispersion relation is exactly 
linear for xa < ^ < 1- The question is thus whether the 
dispersion relation has a local minimum at I — xa, where 
we know the order parameters exactly, i.e., (pAfi — ^1, 
and 4>B,o — ^0- This fact allows us to examine the region 
just below I — XA (and the region just above I — xb, if 
necessary). 

More specifically, if e = — ^ is a small and positive 
quantity, one may argue that Cq cx c| oc e, while dti oa 
d^ (X e. The asymmetry between the two species arises 
because ci — 1 and do — 1 ed I = xa- As a result, for 
component A, cqc\c2 oc c^c^, which implies that C2 oc 
Co, while for component B, d-idgdi oc d^^dQ, and thus 
d-i oc di. All the other coefficients are of higher order 
in e, and thus negligible as / — > x^. Since the stability 
of the persistent currents is determined from the slope 
of the dispersion relation, we may keep only the terms 
which are linear in e. Under these assumptions we find 
that the energy per particle is, up to e, 

Eq/N - 7/2 « / -I- 2xAcl + 2xBd\ + 
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FIG. 1: The critical coupling 7cr given by Eq. (|6]), in the in- 
terval < Z < 1, as a function of xa, for a ring potential. 
The points at = 1 show 7cr for the higher intervals of Z, as 
explained in the text. 



7 [xA{ca + C2) + XB{d-i + di)Y 



(5) 



where we have expressed ci in terms of cq and C2 , and do 
in terms of d_i and di through the normaUzation condi- 
tions. The above expression has to be minimized under 
the constraint of fixed angular momentum, I — xa{c\ + 
2c^)-\-XB{—d'^_i+d\) — XA — i- We do this by minimizing 
the function Eq/N + X[xa{cI + 2c|) + xsi-d'i^ + dj)], 
where A is a Lagrange multipher. The resulting equation 
that connects A, xa, xb and 7 is A(A^ — 4) [A -I- 2{xb — 
Xa)] = 27. For any 7, the above equation has three so- 
lutions, two of which are physically relevant. The one 
appears for < A < 2(xa — xb) — 2{2xa — 1), which is 
< 2, and the other one for A > 2. The first solution gives 
the critical value of 7, 7cr, which gives a zero slope of the 
spectrum Eq/N for < ^ < 1, at Z = as function of 
Xa, namely 



7„. = (3/2)/(4xA-3). 



(6) 



The above expression not only gives the exact value of 
7cr for Xa = 1 and xb — (which is 3/2, as mentioned 
earlier) , but also for any (allowed) value of xa ■ Since the 
above function diverges for xa — > 3/4, persistent currents 
are only possible for 3/4 < a;A < 1. 

In the intervals of higher angular momentum, n < I < 
n+1 with n ^ 0, the situation with stability is rather dif- 
ferent. According to Eq. ([3]) the periodic part of the dis- 
persion relation e{l) repeats itself in each of these inter- 
vals with a slope that is equal to (n-f 1)^ — = 2n+l = 
3,5,7,... For n ^ one has to use the other solution for 
A > 2(xA — Xb)- For the case of only one component, 
Xa = 1 and xb — 0, this solution implies that persistent 
currents are stable for the values 7cr — (2n-|-l)(2n + 3)/2, 
a,t I — n + 1. While the above states support persistent 
currents, as soon as xb becomes nonzero - even if a;^ — > 
but finite - the other solution that lies in the interval 
< A < 2(xa — Xb) has a lower energy, and destabilizes 
the current. In other words, the currents are very fragile 
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FIG. 2: The lowest eigenergies of the Hamiltonian for = 
17, Nb ~ (higher, dashed curve), as well as for Na = 
15, A'^B = 2 (lower, solid curve), with (7 = tt, in the range 
< L < 38, in the truncated space of single-particle states 
with Iml < 7. 



with respect to admixtures of a second species of atoms. 
As a result, the system cannot support persistent cur- 
rents at any interval other than the first one with n ^ 0, 
ior Xb ^ 0. Figure 1 shows 7cr of Eq. ([61), as well as 
the points corresponding to -fci — (2n + l){2n + 3)/2 for 
n = 1,2, and 3. 

To gain some physical insight on the above results, we 
note that for < Z < 1, since the system is in the state 
(pA = <&i and 0s = $0 at Z = XA, it may reduce its 
angular momentum by either transferring some atoms of 
species A from $1 to $0i or some atoms of species B 
from $0 to $-1. However, the second option is energet- 
ically expensive because the angular momentum of $_i 
is opposite to the angular momentum of the system. In 
the second interval 1 < / < 2 (and in any higher one) 
the system is in the state 0^ — ^2 and 4ib = ^1 when 
1 = 1 + Xa- In this case, however, the most efficient way 
for the gas to reduce its angular momentum is to transfer 
atoms of species B from $1 to $0, and not to transfer 
atoms of species A from $2 to $1, as in the first inter- 
val. It is precisely this asymmetry between the first and 
any other interval that allows stable persistent currents 
in the first interval only, but not in any other. 

Beyond the mean-field approximation. To go beyond 
the mean-field approximation, we have also performed 
numerical diagonalization of the Hamiltonian for fixed 
numbers of Na, Nb and L units of angular momen- 
tum. In the case of one component, we have confirmed 
the results derived within the mean-field approximation 
7cr = 3/2 for n = 0, and 7cr = 15/2 for n = 1- 
What is even more interesting is the lowest eigenenergy 
of the Hamiltonian for Na = 17, Nb = 0, as well as 
for Na = 15, Nb = 2, in the range < L < 38, in- 
cuding all the single-particle states with \m\ < 7, for 
U = TT, which is shown in Fig. 2, [the corresponding 
value of 7 has to be calculated according to the formula 
7 = {N—1)U I (27r), which gives 7 = 8]. Figure 2 indicates 
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clearly the metastability of the currents for L = Na, 
and L = 2Na when there is only one component. With 
the addition of even a small second component, the local 
minimum around L — 2Na disappears, destroying the 
metastable current, while the minimum around L — Na 
still exists [7cr ~ 2.83, according to Eq. in agreement 
with the mean-field approximation. 

We have also found numerically that for < L < Nb 
(and Na < L < Na + Nb), the (whole) excitation spec- 
trum is given by the formula Eq{L) = L + ?7/(27r)[q^ -|- 
{N + 1- 2L)q + N{N - l)/2 - L], where g = 0, 1, 2, ... in 
the truncated space of single-particle states with to = 
and 1 (the only ones which are macroscopically occupied 
in the limit of large N). The lowest energy per particle 
Eq{L)/N = I + 7/2 agrees with the result of mean-field 
in the limit N ^ 00, L ^ 00 with L/N = I (finite), and 
NU finite. 

A more specific case of the above spectrum may even 
be derived analytically with use of the Bogoliubov trans- 
formation, for L = Nb (or L — Na), within the same 
truncated space of the single-particle states with to = 
and 1. Within the Bogoliubov approximation, the Hamil- 
tonian takes the form in this case 

H^Nb + U/{27t) [N{N - l)/2+ 

+ iN/2)ia\ai+blbo) + V NANB{aibo + a\bl) ,(7) 

where ai is the annihilation operator of a boson of species 
A with angular momentum to = 1, and 60 is the annihila- 
tion operator of species B with m — 0. This Hamiltonian 
is diagonalized with a Bogoliubov transformation, which 
implies that the eigenvalues are (assuming, for example, 
that Na> Nb), 



£q{NB) = Nb 



U 
2^ 



N 



-{N~2) + {NA-NB){2q + l) 



.(8) 



We then find that the difference Eq{L — Nb) — £q{L — 
Nb) = Uq{q -I- 1) oc and thus vanishes for large N. 

Conclusions. This study provides an interesting illus- 
tration of the physical origin of persistent currents and, 
more generally, of superfluidity. The extra degrees of 
freedom due to the second component, combined with 
the assumed one-dimensionality and the periodicity of 



the Hamiltonian, introduce novel physical effects, which 
have not been known in the physics of the "traditional" 
superffuids. 

More specifically, 

(i) : In one-component systems, sufficiently high values 
of the coupling give rise to persistent currents [l| . In the 
present case, unless the population of the second species 
is sufficiently small - in which case one goes back to the 
one-component case - the second species provides an en- 
ergetically inexpensive way for the system to get rid of 
its circulation: the node that is necessary to form in the 
component that carries the circulation, in order for the 
circulation to escape from the ring, is filled by the second 
component, very much like the coreless vortices studied 
in higher dimensions. 

(ii) The reduced dimensionality introduces another re- 
markable effect: while metastability of persistent cur- 
rents is absent in two-component systems in higher di- 
mensions [1, , here the assumed one-dimensional mo- 
tion makes it possible for persistent currents to be stable, 
at least under specific conditions. 

(iii) The assumed periodicity in the Hamiltonian re- 
flects itself on the dispersion relation, which is quasi- 
periodic, as in the one-component problem. On the other 
hand, while persistent currents corresponding to the flrst 
interval of the angular momentum of the quasi-periodic 
part of the spectrum are stable, for higher values of the 
angular mometum, persistent currents are highly fragile, 
even for a very small admixture of a second species. This 
result is also in sharp contrast to the one-component case. 

The results presented in our study deflnitely deserve 
experimental investigation, in order for our predictions to 
be confirmed. One effect that deserves both theoretical, 
as well as experimental attention is the deviation from 
the one-dimensional motion assumed here. One may ar- 
gue that as this deviation increases, competing mech- 
anisms change the behavior of the system, interpolating 
between one- and two- or three-dimensional motion, thus 
giving rise to rich physical effects. 

Last but not least, in addition to the above more the- 
oretical remarks, the large degree of tunability of the 
persistent currents that we have demonstrated here, also 
makes these systems very appealing in terms of future 
technological applications. 
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